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Abstract 

As a key method in dealing with uncertainties, feedback has been understood 
fairly well in classical control theory. But for quantum control systems, the capability 
of measurement-based feedback control (MFC) has not been investigated systemati- 
cally. In contrast to the control of classical systems where the measurement effect is 
negligible, the quantum measurement will cause a quantum state to collapse, which 
will inevitably introduce additional uncertainties besides the system initial uncertainty. 
Therefore, there is a complicated tradeoff between the uncertainty introduced and the 
information gained by the measurement, and thus a theoretical investigation of the 
capability of MFC is of fundamental importance. In this paper, inspired by both the 
Heisenberg uncertainty principle for quantum systems and the investigation of the 
feedback capability for classical systems, we try to answer the following three basic 
questions: (i) How to choose the measurement channel appropriately? (ii) Is the MFC 
still superior to the open loop control in dealing with the system information uncer- 
tainties? and (iii) What is the maximum capability or limitation of the MFC? These 
questions will be answered theoretically by establishing several theorems concerning 
the asymptotic reachability of eigenstates of a typical class of Hamiltonian control 
mode. 

Index terms Quantum feedback control, quantum measurement, Heisenberg uncertainty prin- 
ciple, impossibility theorem. 
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I. Introduction 

Over the past decades, much progress has been made in regulating the microscopic world, 
such as molecules, atoms and so on (e.g., [1], [2], [3]). The rapid development of the technology 
further promotes the development of the quantum control theory. According to whether or not 
the state information is used in regulating the quantum system during the control process, we 
may divide the control strategy into open-loop control (OLC) and feedback control (FC). For the 
OLC, we can further divide it into Hamiltonian and non-Hamiltonian control modes according 
to how the control law is introduced. For the Hamiltonian control mode (e.g., [4], [5], [6]), we 
mean the control is an adjustable parameter in the Hamiltonian. While in the non-Hamiltonian 
control case (e.g., [7], [8], [9]), one regulates the quantum system by adjusting some parameters 
of an auxiliary system, for example, the parameter of the environment of the system. For the FC, 
there are three types: (i) learning control (e.g., [10], [11], [12]) (ii) measurement-based feedback 
control (MFC) (e.g., [13], [14], [15], [16], [17]) and (iii) coherent feedback control (e.g., [18], 
[21], [22]). For the learning control, one begins with some initial control pulses, and then makes 
some measurement on the samples after the control process. A selected learning algorithm is 
used to adjust the control pulses and then the above steps are repeated on some new samples until 
the performance index does not change significantly. Hence, for this control mode, one needs 
different samples for different cycles. For the MFC mode, one performs some direct or indirect 
measurement on the quantum system and then design a control law based on the measurement 
information to regulate the quantum system. In contrast to the MFC mode, the coherent feedback 
control which was first put forward in [18] does not involve any measurement, and instead, the 
controller is coherently connected with the system plant and can be a quantum system itself. 
This type of feedback has inspired much subsequent research, see, (e.g., [19], [20]). Note that 
what we are really interested in is the state of the system plant. Therefore the controller itself 
will generally cause the quantum decoherence to the controlled system even though it coherently 
entangles with the system plant [21]. Thus whether the effect of the coherent feedback is better 
than the OLC for quantum systems needs to be investigated in depth. 

In this paper, we only focus on the characteristics of the MFC of quantum systems, where 
measurement is needed to gain information of the system to be used for feedback control design. 
For classical systems, the back action effect of the measurement on the system state can be 
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neglected, which, however, cannot be neglected in the quantum case. This is because, by the 
postulates of the quantum mechanics, the measurement on a quantum system will generally cause 
the quantum state collapse, which will inevitably introduce uncertainty to the system besides the 
original uncertainty that the system already shared [23]. Therefore, in the case of the quantum 
MFC, there are some fundamental questions which are not met in the classical control theory, 
for example, 

• How to select the measurement channel appropriately? For example, how to gain the 
information of the system state as much as possible while the corresponding induced 
uncertainty to the system is as little as possible? 

• Even if we have selected a measurement channel appropriately, the measurement induced 
uncertainty is still there, so whether the MFC is still superior to the OLC in dealing with the 
original uncertainty of the system is far from obvious and needs to be investigated seriously 
([24], [25]). 

• Inspired by the well-known Heisenberg uncertainty principle for quantum systems and the 
investigation of the maximum capability of feedback in classical systems (e.g., [26], [27], 
[28]), one may ask what is the maximum capability of the MFC for quantum systems? 

In this paper we will give some partial, yet first rigorous, answers to these three fundamental 
questions respectively. 

The paper is organized as follows. In Section 2, we briefly introduce the typical coherent 
control model and set up the control problems. In Section 3, we give a selection theorem on 
the measurement channel. Then we compare the capability of the OLC and the MFC in dealing 
with the initial state uncertainty in Section 4. We further give an impossibility theorem on the 
control capability of the MFC for quantum systems in Section 5. Section 6 concludes the paper 
with some remarks. 

II. The Coherent Control Model 

We now sketch the control model to be used (see e.g., ([16], [29]) for more details). Let % 
be the Hilbert space of the quantum system, which is of finite dimension dimH = iV < oo. 
The corresponding system state space will be represented by the convex set of density matrices 

S = {pe C N » N :p = p*>0, Tr(p) = 1}, 
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where p* represents Hermitian conjugation of p, C is the complex space of N x N and 
Tr(-) is the trace operator. The set of pure states in S can be represented as 

S = { p e C NxN ■. P = W*, $ e C*, = l}- 

The quantum bath is modeled by the symmetric Fock space T . Now, let A denote the von 
Neumann algebra generated by the set of all bounded operators on HQQJ 7 . Then the quantum 
probability space of the system and the bath is (A,p), where the state p = Pu(£) Pvacuum IS 
given by some state p n on % and the vacuum state p vaC uum on the bath T. 

Under some proper approximations, the quantum dynamics can be described by the Hudson- 
Parthasarathy equation 1 , (e.g., [29], [30]) 

dU t = {-i(H + u t H b )dt - -L*Ldt + yfnLdA* t - ^L*dA t }U t , 

2 (1) 

U = I. 

Here, U t describes the unitary evolution of the whole plant, H , H b and L are system operators, 
and A t is the annihilation operator on T . The symbol * denotes the Hilbert space adjoint as well 
as the scalar complex conjugate. Furthermore, H and H b are Hermitian and can be considered as 
the effective Hamiltonian and the control channel of the system respectively. {u t } is an adjustable 
parameter process and k > is the effective interaction strength. 

In this paper, for the MFC, we will measure the field observable Y t = U*(A t + A* t )Ut- Noting 
the form of the observable Y t and the quantum dynamics (1), L may be called the measurement 
channel. 

For any system observable X, let jt(X) = L^*(X® ldjr)U t . Then from the quantum Ito rule 
and (1) one has 

dY t = y/Kj t (L + L*)dt + dA* t + dA t . (2) 

Now, we want to slightly extend the observation model to contain the case where the detection 
efficiency is not perfect. In fact, there are always some technical noises added to the measured 
signal in practice. We will model these noises by an additional term B t +B% which is uncorrelated 
with A t + A* t and does not interact with the system. Therefore, one has 

dY t = V^VJt(L + L*)dt + y/v(dA* t + dA t ) + \/± - v(dB t + dB* t ), 

'We always work in units such that h = 1. 
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where r\ G (0, 1] is the detection efficiency, with r\ = 1 corresponds to a perfect detection. 
We remark that Y t has the following two crucial properties [16]: 

1) Y t is self-nondemolition, i.e., \Y t , Y s ] = 0, for all s < t. This property implies that the 
observation process can be viewed as a stochastic process on a classical probability space 
by the spectral theorem [31]. 

2) Y t is nondemolition, i.e., \jt(X),Y s ] = 0, for all s < t. This property implies that we can 
estimate the state of the system based on the observations. 

Note that what we are really interested in is the state of the system. Thanks to the above properties 
of the observation process Y t , we can use the quantum filtering theory to estimate the state of 
the system. 

First, we calculate the quantum conditional expectation Tv t (X) = S[j t (X)\B t ] where £ is the 
quantum expectation, B t is the cr-algebra generated by Y s < t . Following the methods of ref. [16], 
we have 

dTr t (X) = 7v t (£X)dt + JrjK(7r t (XL + L*X) 

(3) 

- 7t t (L + L*)ir t (X))(dY t - vV t (L + L*)dt), 
where CX = i[H + u t H b , X] + kL*XL - \{L*LX + XL*L). From the properties of Y t and the 
quantum Ito rule, it can be shown that the innovation process dW t = dY t — v /r/K7r i (L + L*)dt is a 
quantum martingale and dWf = dt [31]. Hence, it is a Wiener process from the spectral theorem 
[31] and the Levy theorem [32]. A rigorous treatment of the quantum stochastic calculus and 
quantum filtering can be found in (e.g., [30], [31], [33]). 

Note that, on the one hand, nt{X) is the conditional expectation of X at time t based on the 
observation of Y s < t , on the other hand, the expectation of X is calculated by Tr[p f X] in the 
Schrodinger picture, so we can define the conditional density operator at time t by 

Tr[p t X] = TTt(X). 

It is easy to get the dynamics of the conditional density operator p t from (3). One has 

dp t = -i[H + utH b , p t ]dt + nV[L]p t dt + ^n[L]p t dW t , (4) 

where W t is a Wiener process on some probability space (fi, J 7 , P) and the superoperators V 
and H are defined by 

V[A}p = ApA* - ^(A*Ap + pA*A), 

DRAFT 



6 



H[A]p = Ap + pA* - Tr(Ap + pA*)p. 

This form of equation (4) is also known as the stochastic master equation or the quantum 
trajectory in physics. In this paper, we do not consider another MFC mode — direct feedback 
control by output measurement, see e.g., [34] and [35]. An interesting comparison between 
indirect feedback control (as in (4)) and the direct feedback control can be found in [36]. 

It is worth pointing out that in the special case where H — 0, H b — F y , L — F z , the equation 
(4) corresponds to the control of quantum spins, which has been investigated in depth, see e.g., 
([16], [17]). While in the current paper, we will focus on more general H , H b , and L to establish 
some fundamental results on the capability of MFC for a wide class of quantum systems. 

To be specific, our control objective is to prepare an arbitrarily desired eigenstate of H 
asymptotically with probability 1 from an arbitrary initial state. Note that preparing a family of 
fiducial states is one of the four basic requirements for quantum computation and is the basis for 
subsequent manipulations, e.g., implementing a family of universal unitary operations on these 
fiducial states [37]. 

Throughout the paper, the admissible control set of MFC is defined as U, which consists of 
control laws adapted to the observation process Y s < t . 

III. A Selection Theorem on the Measurement Channel 

Note that in contrast to the classical control theory, we should choose an appropriate measure- 
ment channel L in the quantum case. This is because the measurement on a quantum state will 
gain some information of the state as well as introduce a state collapse in general. Remember 
that L is called the measurement channel, because it depicts how the back action effect of the 
measurement on the bath is imposed on the system. In this section, we will show how to choose 
the measurement channel L appropriately. 

First, we give the following definition: 

Definition. An eigenstate pf of H is said to be asymptotic reachable under L if there exit 
some control channel H b and control law u t such that for arbitrary initial state p , there exits a 
unique solution of (4) denoted by <^\ L \p G ,H b ,u) 2 , which satisfies 

P{lim <S> ( t L \p ,H b ,u) = p f } = l. 

t— YOO 

2 In the following, it may be abbreviated as &t(po,u). 
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Now, we only consider the case where H and L are non-degenerate 3 , and we have the 
following theorem. 

Theorem 1: Suppose that H and L are non-degenerate, L = L* and [H , L] = 0. Then every 
eigenstate of H is asymptotic reachable under L. 

The proof of the above theorem is based on that of Theorem 4.2 in Ref. [17]. 

First of all, Let us take the N eigenstates ip n , n — 1, . . . , N of H as the basis of the matrix 
representation of the system operators in question, where ip n is the vector with only a nonzero 
element 1 in the n-th row. Then H is a diagonal matrix. Let us define the distance between the 
state p G S and one of the H eigenstates p d = 4>dip*d as 

D(p,p d ) = 1 - Tr {pp d ). 

Secondly, we define 

S >p ±{peS:P<D(p, Pd )<l}, 
S>p±{peS:P<D(p,p d )<l}, 
S<p±{peS:0<D(p 7 p d )<P} 
S <p ±{peS:U<D{p, Pd )<f3}. 

The proof of Theorem 1 will be proceeded by the following five lemmas whose corresponding 
proofs can be found in the Appendices. 

Lemma 1: For the control model (4) with u t being a bounded and real cadlag 4 process and 
p £ S, there exits a unique strong solution p t = Q t (Po, u ) m S. 

Lemma 2: (Proposition 3.5 of Ref. [17]) For the control model (4) with p 6 S and u t = u(p t ) 
where u G C 1 (5, M), there exits a unique strong solution p t = $t(po, w) in S, and u t is bounded. 

3 "Here non-degeneracy of a matrix means that the eigenvalues belonging to different Jordan blocks of its Jordan's normal 
form are mutually different" [38]. In physics, the degeneracy can be lifted by an external field. 
4 Its sample pathes are right continuous and with left limits [39]. 
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Lemma 3: For the control model (4) with u t = 1 , denote A = —i(H + H b ) — nL*L + ay/K,L, 
where a is a real number. If there exist anaeR and an eigenstate pd = ipdipd °f -^o sucn mat 



rank 



N, 



(5) 



in finite 



V r d A N - 1 ) 

then there exists 7 > such that for any initial state p £ <S>i- 7 > p t will exit 5>i_ 7 
time with probability 1. 

Lemma 4: For the feedback control model (4) with u(p t ) = —Tr(i{H b ,p t }pd), denote by 
$t(p, u) the solution of the model with an initial state p and a control u. Then one has 

1-7 



[ sup £>($ t (p,u),A,)>l- 7 /2}< 

0<t<oo 1 



1 — p < 1, for all p E S< 



Lemma 5: For the feedback control model (4) with u(p t ) = —Tr(i{Hb,pt]pd), if H , L are 
non-degenerate, L = L* and [ifo, L] = 0, then the trajectories of p t which never exit <S<i_ 7 /2 
will converge almost surely to pd as t — > 00. 

Proof of Theorem 1. We divide the proof into two steps. In the first step, we show how to 
select the control channel Hb and to construct the corresponding control law u t according to 
the arbitrarily desired eigenstate of H . In the second step, we show that there exists a unique 
solution p t which converges to the desired eigenstate with probability 1 under the designed 
control law u t . 

Step 1. For the feedback control model (4), let us first consider how to select the control 
channel H b appropriately. Let A = —i(H + Hb) — kL 2 + ay/nL as defined in Lemma 3, where 
a is a real number. We proceed to select Hb such that there exists a real number a satisfying 
the following rank condition: 



/ 



rank 



P d A 



\ 



N 



(6) 



for all the eigenstates pd = ipdfl^d (d = 1, - ■ ■ ,N) of H . 

Note that this is similar to the observability condition in the linear control system theory. 
Now, we show how to choose the control channel Hb to meet the above rank condition. 
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First we associate (A; H , H b , L) with a non-oriented graph G(A; H , H b , L) = (V, E), where 
the vertices set V corresponds to the eigenstates of H , i.e., V = {ip d , d = 1, • • • , N}. There is 
an edge between and iff {H b ) i: j ^ for i ^ j. Hence, E = : (Hb)ij ^ for 

i 7^ j} does not depend on L. 

For the feedback control model (4), a necessary condition for A to meet the rank condition (6) 
is that the graph G(A; H , Hb, L) is connected. Actually, Since H is diagonal, non-degenerate 
and [H , L] = 0, we know that L is also diagonal. Thus the off-diagonal elements of A are those 

of —iH b . Hence, if G(A; H , H b , L) is not connected, then we can find a permutation matrix P, 

A l 



such that [40] A = P 



rank 



A 2 



P . Thus for arbitrary d G {1 



/ 



rank 



\ 




Note that i — 1, • • • , N are the basis of the matrix representation, and P is a permutation 
matrix, there exits some integer k 6 {1, ■ ■ • , iV}, which may depend on d and P, such that 

A 2 



rank 



rank 



n 



v 




JV-l 



< N. 



J 



Therefore, the connectivity of graph G(A; H , H b , L) is necessary for A in feedback control 
model (4) to meet the rank condition (6). 

Intuitively, a connected graph G(A; H , H b , L) implies that all the eigenstates of H (the 
corresponding energy levels in physical meaning) can reach each other under the control law u t . 

Moreover, let G\ be the set of the graphs corresponding to matrices H b whose elements on 
the secondary diagonal are all nonzero while the others are all zero. From the proof of Lemma 
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4.4 in [17], it is easy to see that if G(A; H , H b , L) G G\, then there exists a such that A meets 
the rank condition (6). Therefore, for the feedback control model (4), a sufficient condition for 
A to meet the rank condition (6) is to choose the control channel Hb such that G(A; H , Hb, L) 
is actually a path 5 of length N — 1. 

Next, we give the feedback control law u t . 

In order to globally prepare an arbitrary eigenstate p d of H with probability 1, let 7 > be 
as defined in Lemma 3, and define B = {p : 7/2 < Tx(pp d ) < 7}, we construct the feedback 
control law as follows: 

1) If Tx(p t p d ) > 7, ut = -Tr(i[H b , p t ]pd); 

2) If Tt{p tPd ) < 7/2, u t = 1; 

3) If p t G B, then u t = —Tr(i[H b , pt]pd) if Pt enters B through the boundary Tv(p t p d ) = 7; 
Ut — 1 otherwise. 

Step 2. In this step, we combine the results from Lemma 1 to Lemma 5 to show the existence, 
uniqueness and global stability of the solution of the model (4) under the control law u t . 

As defined in the proof of Lemma 3, $t(p, u) denotes the solution of equation (4) at time t 
with the initial state p and the control u. Now, for an arbitrary fixed initial state p G S and the 
control law u given in step 1, we first construct a solution $ 4A „(p, u) , where n is a positive 
integer. 

Define the stopping time 

7? = n A n = mf{t > : $ t (p, 1) G <S< X _ 7 } A n. 

Then $ tAn (p, 1) = $ f (p, 1), for t < r" and p T n = $ T ™(p, 1). In the following, we denote 
$s,t(Ps 5 M ) as me solution of equation (4) at time t with the initial state p s at time s and the 
control it, and define the stopping time 

(r"=(TiAn = inf {t > rf : $ r ™ t (p T », u) G <S>i_ 7/2 } A n, 

where u = — Tr(z[if&, Pt\Pd)- Then we know for t < cr", 

$tAn(p, «) = ^<r™^(p, 1) + ^r<t«rr^Tr,t(pTr,«), 

5 A path in G(A; Ho, Hi,, L) of length r is a sequence [ipi , if>i r ] of distinct vertices such that (tpi i _ 1 , tpij ) € E. 
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where Xa represents an indicator function on the set A. Then we can recursively define a 
sequence of stopping times as follows 

ol = a k A n = mf{t > r k : $ T g, t (p T n,u) G <S>i_ 7 / 2 } A n, 

and 

T nA TkAn A inf { t > . $^^0^, 1) G S<^} A n, 
where p a n = $ r n j(J n(p T n, u) and p T n = $<t%_ 1 ,t , q{.Po%_ 1 i !)• Now we can construct the solution 

oo 

®tA*(j>, U ) = Xt<T?®t(P, 1) + ^2( X t£ <t«7» ,t(Pr^ , «) + ^<t<7£ +1 $<x£,t(Ax£, 1)) 

fc=l 

for t<a n = crAn = linifc^oo a k An < n. Since from Lemmas 1 and 2, every segment between 
any two of the stopping times is a.s. uniquely defined, it is easy to know that the solution is 
a.s. unique on t G [0,cr rt ]. Now let n — > oo, we obtain the unique solution $ t (p, w) defined up 
to time a. Furthermore, we need to prove that the solution exists for all t > 0, which can be 
derived just by showing a = oo a.s.. Below we will show that only finitely many a k are finite 
for almost every sample path. 

We can get the fact that the strong Markov property holds on each segment between any two 
of the stopping times r k < t < o k and a k < t < T k +\ from Proposition 3.7 in [17]. Hence by 
Lemma 4 we know that 

P{(Tfc < Oo|Tfc < OO} < 1 — p, 

and by Lemma 3 we have 

¥{r k < oo|cr fc _ 1 < oo} = 1. 

Moreover, 

F{r k < oo|crjt < oo} = P{o"fc_i < oo\rk < oo} = 1. 

Hence we have 

P{<7 fc < oo} _ F{r k < oo\a k < oo}P{cr fc < oo} P{crfc_i < oo|r fc < oo}P{r fc < oo} 
PK_x < oo} ~ P{r fc < oo} P{a fe _! < oo} 

= P{cr fc < oo|rfc < oo}P(r fc < ooluk^x < oo) 
< 1 -p. 
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Hence, we have 



F{a k < oo} < (1 -p) 



k-l 



SO 



oc 



00 



1/p < 00. 



k=l k=l 

Therefore, by the Borel-Cantelli lemma [41], we know that 



P{ lim <7fc = cr = 00} = 1. 



Thus for almost every trajectory, there exists a positive integer K{u) < 00 such that = 00 
for any k > K, and 07. < 00 for any k < K. Then from Lemma 5, we have that $t(p, u) — > pd 



In this section we will compare the effect of the quantum OLC and MFC in dealing with 
the initial state uncertainty. It is well known that in the classical control theory, the feedback 
control is much superior to the OLC in dealing with uncertainties, such as parameter uncertainty, 
structure uncertainty, external disturbance, or all, see e.g., ([26], [27]). However, the comparison 
of the quantum OLC and MFC is far from obvious. This is because, as we have stressed before, in 
the measurement-based quantum feedback control, the measurement generally introduces another 
kind of uncertainty on the system besides its initial uncertainty. Hence, during the control process, 
we have to deal with these two kinds of uncertainties. Actually, we have proposed and analyzed 
this question in [24] based on a special model, and demonstrated that the measurement-based 
quantum feedback control is still superior to the OLC in some sense. 

The measurement-based quantum feedback control model is represented as equation (4). As 
we have mentioned, when referring to the OLC strategy, we only use the prior information to 
design the control law. There are two types of OLC models: (1) there is no measurement at all; 
(2) there is measurement but is not used to design the control law. 

Our control target is also to prepare an arbitrarily desired eigenstate of H from an arbitrary 
initial state. 

First, we look at the case where there is no measurement. The OLC model is described by 



a.s., as t — > 00. The proof of Theorem 1 is completed. 



□ 



IV. Comparison with Open-loop Control 



dpt 
dt 



i[H^+u(t)F yiPt } 



(7) 
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where H' may be different from the effective Hamiltonian H Q in equation (4) due to the 
measurement back action effect originally introduced there [29]. 

Theorem 2: For the OLC model (7), for arbitrary control channel H b and arbitrary control 
law u t , an arbitrary eigenstate of H cannot be prepared from any mixed initial state. 
Proof. It is easy to see that the evolution of equation (7) is unitary. Consider the von Neumann 
entropy 

S(p) = -Trfplogp], 

which has the following two properties: 

1) It is invariant under the unitary evolution; 

2) The entropy is zero if and only if the state is pure. 

Therefore, if the initial state is mixed, its entropy S(p ) > 0, and then S(p t ) = S(p ) > for 
all time t. Remember that the target state is an eigenstate of H , which means its entropy is 0. 
Hence, we cannot prepare the target state from a mixed initial state no matter how to choose 
the control channel H b and how to design the control law u t by using the OLC model (7). □ 

Note that the initial state of the system may generally be a mixed state 6 because of the 
inevitable interactions between the system and the bath. Hence, we cannot achieve the target no 
matter how to choose the control channel H b and how to design the corresponding control law 
u t if we use the OLC model (7). 

Next, let us look at the case where there is measurement but only the prior information is 
used to design the control law. The corresponding OLC model is [29] 

^ = -i[H + u t H b , pt] + KV[L] Pt . (8) 

In contrast to equation (4), there is no diffusion part in equation (8) which depicts the mea- 
surement induced uncertainty. This is because equation (8) actually describes the evolution of 
the ensemble average. This form of evolution equation is usually called the master equation in 
physics. We have the following two theorems. 

Theorem 3: For the OLC model (8), if H is non-degenerate, and [H , L] = 0, then for 
arbitrary control channel H b and arbitrary control law u t , one cannot prepare an arbitrarily 
desired eigenstate of H from any mixed initial state. 

6 The von Neumann entropy of a mixed state is strictly greater than 0. In this sense, we say the initial state has some 
uncertainties. 
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Proof. Let us first consider the dynamics of Tr(p 2 ). From equation (8), we have 

dTrjp 2 ) 
dt 

= 2Tr[-i(H + u t H b )p 2 t + tp 2 (H + u t H b ) 

+ KLp t L*p t -^(L*Lp 2 t +p 2 t L*L)] 
= 2nTr[Lp t L*p t - L*Lp 2 ] 

< 2KT^[Lp 2 L*]T^[LL*p 2 } - 2KTr[L*Lp 2 } 

< K(Tr[Lp 2 L*} + Tr[LL*p 2 ]) - 2kYx[L* Lp 2 ] 
= Klr[(LL* - L*L)p 2 }. 

Since H is Hermitian and non-degenerate, it is a diagonal matrix and its eigenvalues are mutually 

different. From this and [H , L] — 0, L is also a diagonal matrix, and so LL* — L*L = 0. Hence, 

^Tr( Pf 2 ) 
dt 

which implies that for all t > 0, Tr(p 2 ) < Tr(pg). Note that Tr(p 2 ) = 1 iff p is a pure state. 
Hence, we conclude that the desired target cannot be prepared from a mixed initial state p 
under the conditions of the theorem. □ 
Theorem 4: For the OLC model (8), if H is non-degenerate, and [H , L] ^ 0, then there is at 
least one eigenstate of H denoted as p d , such that no matter how to select the control channel 
H b and how to design the corresponding control law u t , one has 

lim sup D(p t ,p d ) >S d >0, 

t— S-OO 

where 5 d = Ti 2 (p d V[L]p d )/{2[2T^(L*L) 2 + Tr^(L*p,L) 2 + Tr^ (L* LL* Lp d )} 2 } . 

Proof. Since H is diagonal, non-degenerate and [H ,L] ^ 0, we know that there exists at 
least one non-diagonal entry of L which is not zero. Without loss of generality, suppose it is in 
the d-ih. column of L. Hence we have 

N 

Ti( Pd V[L] Pd ) = -Y,L; d L id <0. (9) 

Next we use a contradiction argument to give the proof of this theorem. Suppose that for any 
eigenstate pi of H , their exist a corresponding control channel H b and a corresponding control 
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law u t , such that 

limsupD(p t ,pi) < Si. 

t— >oo 

Then there exist some e > 0, T > 0, such that for any t > T 

D(pt,Pi) <5i-e. 

Therefore for t > T, 

\\pt-Pi\\ =TrH(pt-p l ) 2 ) < (2D(p uPi )^ < y/2(6i-e). (10) 
Then by (10) we have 

K(Ti(p t V[L]p t ) - Ti(p d V[L]p d )) 
= K(Tr(Lp t L* Pt ) - Tr(Lp d L*p d ) + Tr(L*Lp 2 ) - Tr(L*Lp 2 )) 

= hi(Tr(Lp t L*(p t - p d )) + Tr((p t - p d )L*p d L) + Tr(L*Lp d (p d - p t )) + Tr(p t L*L(p d - p t ))) 

< K(Tr^(Lp t L*f + Tt^ (L* p d L) 2 + Tr^ (L*LL*Lp d ) + Tr^ (L*LL*Lp 2 ))\\p t - p d \\ 

< ^/2(5 d -e)K(2T^(L*L) 2 + Tr^ (L* p d L) 2 + Tr^ (L*LL*Lp d )). 

(ID 

Then by (9) and (11), we have for t > T 

^1 = 2^{ Pt V[L]p t ) 

= 2KTr( Pd V[L} Pd ) + 2K(Ti(p t V[L]p t ) - Ti(p d V[L]p d )) 

< 2KTi(p d V[L] Pd ) + 2^2(5 d -e)K(2Ti l *(L*L) 2 + Tr^ (L*p d L) 2 + T^(L*LL*Lp d )) 
= a(e). 

By substituting the definition of 5 d into a(e) defined above we can find that a(e) < 0. Hence, 
from the above inequality, we have a contradiction 

< Tr(p^) -»■ -oo, as t -> oo, 

and the proof is completed. □ 
From Theorems 2, 3 and Theorem 4, we know that if using the OLC model, we cannot achieve 
an arbitrarily desired target state no matter how to select the measurement channel L, the control 
channel Hf, and the control law u t . Hence, in comparison with the MFC results as established in 
Theorem 1, we conclude that the measurement-based quantum feedback control is still superior 
to the quantum OLC in dealing with the initial state uncertainty. 
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Note that in contrast to the classical control theory, the measurement-based quantum feedback 
control model (4) is different from the OLC models (7) and (8). This is due to the inherent 
quantum measurement back action effect which consists of the deterministic drift part and the 
uncertainty part. Specifically, the change from H' in (7) to H in (4) and nV[L]p t dt describe the 
deterministic back action effect, while ^JWjl-L[L\p t dW t depicts the uncertainty part. It is these 
structural changes, especially the uncertainty part, that make it possible to design a feedback 
control law to achieve the control target. 

V. The Impossibility Theorem of MFC 

In this section, we will give an impossibility theorem on the MFC if the selected measurement 
channel L and the effective Hamiltonian H do not commute. 

Theorem 5: For the MFC model (4), if H is non-degenerate, [H ,L] ^ and 77 e [0,1), 
then there is at least one eigenstate of H denoted as pd, such that no matter how to select the 
control channel H b and how to design the corresponding control law u t , one has 

lim sup .D^pd) > A d > 0, a.s. 

where A d = [2Tr(p d V[L]p d ) + V Ti(n[L]p d ) 2 ] 2 /{2[2cp 1 (L, p d ) + v <p 2 (L, Pd )] 2 }, 
<Pi (L, Pd ) = 2Tr I (L*L) 2 + Tri (L*p d L) 2 + Tr I (L*LL*Lp d ) , 
<p 2 (L lPd ) = 2Tr^(L*L) 2 + 2Tt^ (L*LL*Lp d ) + 2Tr(L*L) + 2Tr^ (LL*p d L*Lp d ) 
+ 3Tr^(L + L*) 2 Tr((L + L*)p d ) + 3Tr^(L + L*) 2 Ti^(L + L*) 2 
+ 2Tr^(L + L*) 2 Tr5((L + L*) 2 p d ). 
Proof. The proof of this theorem is similar to the proof of Theorem 4 and a contradiction 
argument will be used. First, since H is diagonal, non-degenerate and [H , L] ^ 0, we know 
that L has at least one non-zero non-diagonal entry. For simplicity, we suppose it is in the d-th 
column of L. Thus for the eigenstate p d , we have 

N 

2 K Tr(p d V[L]p d ) + r] K TT(n[L]p d ) 2 = -2«(1 - 77) ^ KdUd = -2«(1 - r))a < 0. (12) 

Suppose that for any eigenstate pi of H , there exist corresponding control channel Hb and 
control law u t , such that 

limsupD(p i7 pi) < Ai, on A, 

t— >oo 
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where A is a set with positive probability F{A} > 0. Then for each sample point u e A, there 
exist e > and T{uS) > 0, such that for any t > T(u) we have 

D(p t (u)),pi) < Ai - e. 

Therefore, for t > T(u) with cu G A, we have 

\\pt{co) -Pif < 2D( Pt (u),pi) < 2(Ai-e). (13) 

Then for each cu G A, by (13) we have 

K(Tr(p t V[L]p t ) -Tr(p d V[L]p d )) 
= K(Tr(L Pt L*p t ) - Tr(Lp d L*p d ) + Tr(L*Lp 2 ) - Tr(L*Lp 2 )) 

= hi(Tr(Lp t L*(p t - Pi )) + Tr((p t - p d )L*p d L) + Tr(L*Lp d (p d - p t )) + Tr(p t L*L(p d - p t ))) 

< k(Tt^ (Lp t L*) 2 + Tt^ (L* Pd L) 2 + Tr^(L*LL*Lp d ) + T^(L* LL* Lp 2 t ))\\Pt ~ Pd\\ 

< ^2(A d - e) K (2Tr? (L*L) 2 + T^(L*p d L) 2 + Tr^ (L*LL*Lp d )) 
= Kipi (L, p d ) y / 2(A d - e), 

(14) 

and similarly for u G A 

TiCH[L]p t ) 2 -TTCH[L]p d ) 2 
= 2Tr(L*Lp 2 t ) - 2Tr(L*Lp 2 ) + Tr(Lp t Lp t ) - Ti(Lp d L Pd ) + Tr(L*p t L*p t ) - Ti(L*p d L*p d ) 
- 2Tr((L + L*) Pt )Tr((L + L*)p 2 ) + 2Tr((L + L*) Pd )Tr((L + L*)p 2 ) 
+ Tr 2 ((L + L*)p,)Tr(p 2 ) - Tr 2 ((L + L*)p d )Tr(p 2 ) 

< 2(Tr^(L*L) 2 + T^{L*LL*L Pd ))\\ Pt - p d \\ + Tr(Lp t L(p t - p d )) + Tr((p t - p d )Lp d L) 
+ Tr{L*p t L*{p t - Pd )) + Tr((p t - p d )L*p d L*) + 2Tr((L + L*)p d )Tr((L + L*)(p d - p t )) 
+ 2Tr((L + L*) Pt )Ti((L + L*)p d (p d - p t )) + 2Tr((L + L*)p t )Tr((L + L*)(p d - p t )p t ) 
+ Tr((L + L*)(p d + p t ))Tr((L + L*)(p t - p d )) 

< {2Tr^(L*L) 2 + 2Tr^ (L*LL*Lp d ) + 2Tr(L*L) + 2Tr^ (LL*p d L*Lp d ) 
+ 3Tr^(L + L*) 2 Tr((L + L*)p d ) + 3Tr^(L + L*) 2 Tr^(L + L*) 2 

+ 2Tr5(L + L*) 2 Tr^((L + L*) 2 p d )}\\ Pt - p d \\ 
= (f2(L, Pd )\\p t - p d \\ < (p 2 (L,p d ) ^2(A d - e). 

(15) 
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By the It 6 formula, we have 

dTr(p 2 ) = 2KTi( Pt V[L] Pt )dt + V KTi(H[L]p t ) 2 dt + 2^/rjKTr{p t H[L] Pt )dW t 

= (2KTi(p d V[L]p d ) + ^Tr(H[L]p d ) 2 )dt + (2KTr(p t T>[L]p t ) + V KTT(H[L]p t ) 2 
- 2nTi(p d V[L]p d ) - r)KTr(H[L]p d ) 2 )dt + 2^rfKli(p t H\L\p t )dW t . 

(16) 

Then by (12), (14), (15) and (16), for t > T(u) with u e A, we have 

Tr(p t 2 ) = Tr(p 2 () ) + £ 2 K Tr(p d P[L] P(i ) + r^Tr(H[L]p d ) 2 ^ 

+ j£ 2 K Tr(p s P[L]p s ) - 2«:Tr(p d P[L]p d ) - V KTi(H[L}p d ) 2 + p K Tr(^[L]p s ) 2 d S 
+ 2VP/ t *Tr(p s ^[L]p s )« 
< Tr(p 2 J + {-2k(1 - n)a + 2k^(L, p d )v/2(A d -e) 

+ p K (p 2 (L,p d )v/2(A d -e)}(t - t ) + 2^ /* Tr ( Ps H[L] Ps )dW s . 



(17) 



Moreover, by Theorem 1.51 in [42], for < e < |, we have 



Tr(p 2 ) < Tr(p 2 J + {-2«(1 - p)a + 2k ¥ > 1 (L, p^v^A.-e) 

+ pK<p 2 (L, p d )v/2(A d -e)}(t - t ) + o{U^) + 0(1). 
— > — oo, as £ — > oo, on A 
where the limit is derived by the definition of A^. This is a contradiction. Thus the theorem is 
proved. □ 
Note that if [H , L] = 0, then all the off-diagonal elements of L must be zero, so by eq. (12), 
we have A rf = 0. Hence, the above limit to the MFC is due to the non-perfect detection and 
the noncommutative relationship between the effective Hamiltonian H and the measurement 
channel L. Theorem 5 is actually inspired by the Heisenberg uncertainty principle which says 
that given two observable H and L, for any quantum state p, their standard deviations satisfy 

u _ 2 

where the mean value of H is (H ) = Tr[H p] an d the standard deviation of H is AiJ = 



\/ (Hq) — (H ) 2 . We see that the non-commutability plays an important role in the Heisenberg 
uncertainty principle, just as the role played in the asymptotic reachability of Theorem 5. 
Combining Theorem 1 and Theorem 5, one has 

Theorem 6: Suppose H and L are non-degenerate, L = L*, and p E (0,1). Then every 
eigenstate of H is asymptotic reachable under L if and only if [H , L] — 0. 
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VI. Concluding Remarks 



In this paper, three fundamental problems concerning the asymptotic reachability of eigenstates 
by the MFC are investigated. First, we have given a selection theorem on how to select the 
measurement channel. Secondly, we have compared the MFC and the OLC in dealing with the 
system initial uncertainty, and have shown that the MFC is still superior to the OLC for the 
quantum systems in spite of the additional uncertainty induced by the measurement. Thirdly, we 
have given an impossibility theorem on the capability of the MFC. It is worth emphasizing that 
the study on the capability of MFC is just initiated. For future investigations, it may be necessary 
to further study the limitations of quantum MFC, e.g., in dealing with structural uncertainties 
and/or in achieving other control objectives. Also, it would be interesting to study how the 
measurement channel and/or the control channel can be adjusted adaptively according to the 
quantum state in question, as discussed in e.g., [44]. The study of the above problems with the 
non-Markovian model would also be interesting. There is no doubt that these investigations will 
help us understand more about the MFC. 

VII. Appendices 

In this appendix, we will give the proofs of Lemmas 1, 3, 4 and 5. The analyzes are mainly 
based on those in [17] but for the more general model in the paper. 

A. Support Theorem 

The following support theorem from e.g. ([17], [39]) connects the properties of solutions of 
a stochastic differential equation to solutions of an associated deterministic system. 

Theorem 7: Let M be a connected, paracompact C°° -manifold and let Xk, k = 0, ...,n, be 
C°° vector fields on M such that all linear sums of X k are complete. Let X k = X l k (x)di in 
local coordinates and consider the Stratonovich equation 



n 




k=l 



Consider in addition the associated deterministic control system 



d 



n 



< = x (<) + ]Tx fc (<)o M fc (t), 




k=l 
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with u k G the set of all piecewise constant functions from 1Z + to 1Z. Then 

y x = {x u : u i c 

where ^ is the set of all continuous paths from 1Z + to M starting at x, equipped with the 
topology of uniform convergence on compact sets, and 5? x is the smallest closed subset of W x 
such that F{lu G Q, : x.(uj) G ^} = 1. 

B. Proof of Lemma 1 

Proof. Consider the corresponding unnormalized linear form of equation (4) 

dp t = -i[H + u t H b , pt]dt + n(Lp t L* - ^(L*Lp t + p t LL*))dt + y/^(Lp t + p t L*)dY t . (18) 

By [39] (pp. 195- 197), it is not difficult to check that the above equation has a unique strong 
solution as it obeys a functional Lipschitz condition. 

Now we introduce a process W t = y/rftt + y/T^rjW, where W t is a Wiener process which 
is independent with Y t . Then with the same initial condition p = po(= Po)> the process p t is 
stochastically equivalent to E[p t |J^], where p t satisfies the following equation 

dp t = -i[H + u t H b , p t ]dt + n(Lp t L* - l -{L*Lp t + p t LL*))dt + yfk{Lp t + p t L*)dW t . (19) 

If p = ipoipo* is a pure state, it is easy to check by the Ito formula that the solution of the 
equation 

dip t = ~i{H + u t H b )$ t dt - ^L*L?p t dt + ^Lip t dW t (20) 

satisfies p t = ^t"4>t • From [39], pp. 326 we know that tjj t ^ a.s. for tjj ^ 0. Therefore, 
p~t = ^\.Pt\Ft \ > a.s.. Note that equation (18) and (19) are both linear, hence, if the initial 
state po is mixed, i.e., p = J2i ^i^oV'o* w i m me convex weights A, > 0, we can also get 
p~t = | Tj\ > a.s.. This implies that Trp t > a.s.. Now it is easy to prove that equation 
(4) with u t G U has a unique solution p t = p t /Tr(p t ) G S by the Ito formula. □ 

C. Proof of Lemma 3 

Proof. First of all, from Lemma 1, we know that equation (4) with u t = 1 has a unique 
solution in the compact set S. Now for some fixed p d = ipcfll>d> we define the function 

g(p) = min E D(<f> t (p, l),p d ), p G S, 
te[o,T] 
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where T is a fixed time and $t(p, 1) denotes the solution of equation (4) at time t with initial 
state po = p and control u t = 1. The proof of Lemma 3 is divided into three steps. 

Step 1. We first show that g(p) < 1 for any p G Si = {p G S : D(p, p d ) = 1}. 

A contradiction argument is now used to complete the proof. Suppose that g(p) = 1 for some 
p G Si, i.e., ETr($ t (p, l)p d ) = 0, for some p G Si and all t G [0,T]. 

From the proof of Lemma 1, we can write $ t (p, 1) = p t /Tr(p t ) with p t = ^\ AjE^^lJ^], 
where Aj's are the coefficients of the spectral decomposition of p = p = J^i-^'oV'o*' e 
T^/^lO} and tjj l t obey the equation 

#j = -i(fio + #b)$d* - 7:L*Li/jidt + ^LtPidW tj (21) 

Now, one can show that ETr($ t (p, l)p<i) = is equivalent to 

ETr(p> d ) = ^Ef^/)*^/)] = 0, 

i 

i.e., g{p) = 1 iff = 0} = 1, for all i and t G [0, T]. 

Remember that 

dW t = ^jdY t + y/T^jdW t 

= r/v^Tr[(L + L*)p t ]dt + 07<fW t + ^1 - rjdW u 

where W t and W t are independent Wiener processes. Since Ti[(L + L*)p t ] is bounded, by 
Girsanov Theorem [32], there exists a measure Q, such that the stochastic process W t is a 
Wiener process in [0, T] under the new measure Q and Q is equivalent to the measure P. Hence, 

g(p) = 1 iff ®{i(>$t = 0} = 1, for all i and t G [0, T]. (22) 
Moreover, it is easy to find that the linear Stratonovich form of equation (21) is 

di)\ = -i(H + H b )4>\dt - KL*Li)\dt + ^JkL^X o dW t , (23) 

i>l = 

Now we apply Theorem 7. Recall that denotes the set of all continuous paths starting at t/Sq 
and yjj o denotes the smallest closed subset of satisfying that Q{cu G f2 : $ l t (uS) G ^L} = 1- 
Moreover, we denote ( = {v t G #^- o : f *^ = 0}. It is easy to prove that f is closed in 
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the compact uniform topology for any t G [0, T\. Then it is obvious that C 5?% t for all i 
and t G [0,T] from (22). 

Now consider the following deterministic differential equation 



(it 



A$ = (-i(ff + #&) - + aVKL)if>l = 



(24) 



Then by Theorem 7 we know that the solution of equation (24) with some a satisfying 
equation (5) belongs to =5^, which implies ip\ G »^ 0it , for all t G [0, T]. That is to say, 

^\*^ d = o, for all t G [0, T]. (25) 

However, from equation (24) and (25) we have 

d 



— WMt=o = (r d A n )% = 0, n — 1,2, ... N 



hence, 



PA 



1>l 



o. 



V PA"' 1 J 

Consequently, from (5) we know that = 0, a. contradiction, and the proof of step 1 is 
completed. 

Step 2. We show that there exists 7 > such that g(p) < 1 — 7 for all p G <S>i_ 7 . 

We first prove that there exists x > such that g(p) < 1 — x for all p G 5>i_ x -. Otherwise, 
suppose that for any x > there exists p x G iS>i_ x such that 1 — x < g(p x ) < 1. 

Take a sequence x n I such that the corresponding sequence p Xn G 5>i_ a;n and g(p Xn ) 1. 
Because of the compactness of 5, we can take a subsequence denoted also as x n such that 
Px n — > Poo £ From step 1 we have g(poo) < 1> i-e., there exist s G [0, T] and e > such that 

EL>($ s ( Poo ,l),p d ) = 1-e. 

Because the solution of equation (4) is Feller continuous by Proposition 3.6 in [17], we have 

1= lim g(p Xn )< \imED(<l> s (p Xn ,l),p d )=ED(<f> s (p 00 ,l),p d ) = l-e<l, 

n— >oo n—>oo 

which is a contraction. Hence, there exists x > such that g(p) < 1—x for all p G <S>i_a;. Take 
7 = x/2, then the proof of step 2 is completed. 
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Step 3. In this step, we will give the proof of this lemma. Let r p be the first exit time of 
1) from <S>i_ 7 with the initial state p £ <S>i_ 7 . To complete the proof of this lemma, we 
just need to show that 

sup Et p < oo. (26) 

p£S>i-7 

Actually, from Lemma 4.3 in [43] we have 

< T 
Tp ~ l-sup e65 P{r,>T}- 

Then in order to prove (26), we just need to show that 

supP{r e > T} < 1. 

Obviously the inequality holds if g £ «S\«S>i_ 7 . Let us just consider g £ <S>i_ 7 . We use a 
contradiction argument. Suppose that for every e > 0, there exists g(e) £ <S>i_ 7 , such that 
P{T e ( e ) > T} > 1 — e. Then we can get 

E£>($ s (f?(e), 1), p d ) > (1 - 7 )(1 - e), for all s £ [0, T]. 

Now we take a sequence e n 4- 0, then we can find a corresponding subsequence denoted also as 
g(e n ) such that g(e n ) 4- £?(oo) £ 5>i_ 7 as n — ^ oo by the compactness of S. Hence by the Feller 
continuity of the solution of equation (4), we have 

ED($Moo),l),Pd) > (1-7), for all s £ [0,T], 

which is in contradiction with the result of Step 2. Hence the inequality (26) holds, and the 
proof of this lemma is completed. □ 

D. Proof of Lemma 4 

Proof. From Lemma 2, we know that equation (4) with u(p t ) = — Tr(i[H , pt]pd) has a unique 
solution in the compact set S. By the Ito formula, we can get 

dD(jH, p d ) = -u 2 t dt - y/rfKTi[(Lp t + p t L* - Tr[(L + L*)p t ]p t )p d ]dW t . (27) 

Denote C as the the weak infinitesimal operator of $<(p, u), then from (27) we have 

CD(%(p,u),p d ) = -u 2 (p)<0, 

from which the proof can be completed by Theorem 2.2 in [17]. □ 

DRAFT 



24 



E. Proof of Lemma 5 

Proof. First of all, from Lemma 2, we know that equation (4) with u(p t ) = — Tr(i[H , pt]pd) 
has a unique solution in the compact set S. We divide the proof into two steps. In the first step, 
we will show that the trajectories of p t which never exit iS<i_ 7 / 2 converge in probability to p d . 
Then we will complete the proof of this lemma in the second step. 

Step 1. Consider the following Lyapunov function 

V(p) = l-Tr 2 (p Pd ). 

It is easy to see that V(p) > in S and V(p) — iff p — p d . 

Note that H is diagonal, non-degenerate and [H , L] = 0, then we know that L is also 
diagonal. In addition, L is non-degenerate, so we know that L and H have the same eigenstates. 
Then from the ltd formula, 

dV(p t ) = (-2u 2 Tr(p t p d ) - A V K(\ d - Tr(L Pt )) 2 Tr 2 (pp d ))dt 

(Z5) 

- 2 v /^Tr((Lp t + p t L* - Tr[(L + L*)p t ]p t )p d )Ti(p t p d )dW t . 
Denote by C the weak infinitesimal operator of p t , then from (28) we have 

CV(p) = -2u(p) 2 Ti(pp d ) - A V K(X d - Tr(Lp)) 2 Tr 2 (pp d ) < 0, (29) 

where is the eigenvalue of L corresponding to the eigenstate p d . Hence from Proposition 
3.4 and 3.7 in [17] and the stochastic stability theory (Theorem 2.3 in [17]), we know that p t 
converges in probability to the largest invariant set contained in J\f s = {p E S : CV(p) = 
0}. Since we just consider the trajectories that never exit <S><i_ 7 / 2 , from (29) we know that 
Afs H 5<i_ 7 / 2 must be a subset of the largest invariant set M = {p E S : Tr(Lp) = A^}. 
Moreover, the invariant set M can only contain p such that Tr($ t (p, u)L) is constant. Therefore, 
by the Ito formula, 

= dTr(Lp t ) = -iu t Tr([H b , p t ]L)dt + 2 v ^(Tr(L 2 p < ) - Tr 2 (Lp t ))dW t . 

Hence we must have 

Tr(L 2 p) - Tr 2 (Lp) = 0. 

Since L is non-degenerate, we conclude that p must be the eigenstate of L. Therefore, it is easy 
to see that As H 5<i_ 7 / 2 = {p d }, so we know that the trajectories that never exit S<i_ 7/ / 2 must 
converge in probability to p d . 
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Step 2. In this step, we will give the proof of this lemma. Let us define the even B = {u G 
Vt : p t never exits <S<i- 7 /2}- Then from step 1 we have 

lim F{D{p t , p d ) > e\B} = 0, for all e > 0. 

t— >oo 

Since D(p, p d ) < 1, for all p £ S, we know that for all e > 

E[D(p t ,p d )\B] <¥{D(p t ,p d ) > e\B} + e(l-¥{D(p t ,p d ) > e\B}). 

Hence, 

limsupE[D(p t ,p d )\B] < e for all e > 0. 

t— >oo 

Thus, it is easy to see that 

UmE[D(p t ,p d )\B] = 0. 

t— >oo 

Note that CD(p,p d ) = —Tr 2 (i[Hb, p]p d ) < in 5<i_ 7 /2, then by the Theorem 2.2 in [17], 
D(p t , p d ) converges a.s. for path remaining in 5 < i_ 7 / 2 . Therefore, by the dominated convergence 
theorem we have 

E[lim D{p t , p d )\B] = 0. 

t— >oo 

So 

P{lim D(p t ,p d ) =0\B} = 1, 

t—>-oo 

which completes the proof of this lemma. □ 
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